A Model for Striped Growth 
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We introduce a model for describing the defected growth of striped patterns. This model, while 
roughly related to the Swift-Hohenberg model, generates a quite different mixture of defects during 
phase ordering. We find two characteristic lengths in the system: the scaling length L(t), and the 
average width of the domain walls. The growth law exponent is larger than the value of 1/2 found 
in typical point defect systems. 



I. INTRODUCTION 



Our understanding of the growth of striped patterns 
after a temperature quench from an isotropic state re- 
mains limited [1-8] . We do not have the simple standard 
scenario of point defect annihilation observed in systems 
like the XY model [9] where there is ordering ending in 
a homogeneous phase. Instead we have pattern coarsen- 
ing via a competition between point and line defects con- 
structed from building blocks of dislocations and disclina- 
tions. Experiments and numerical simulations for simple 
models are in agreement that there is scaling in stripe 
forming systems governed by a characteristic length or 
growth law, L(t) w t x , with an exponent x « 1/4 to 1/3 
[8] . This exponent is smaller than one would expect from 
the simplest theoretical treatments which give x ~ 1/2. 
A convincing theoretical understanding of the values of 
x for isotropic pattern forming systems is still lacking. 
A complication is that there may not be a well defined 
value for x, but, as found in numerical treatments, x may 
depend weakly on a control parameter. 

When we turn to the defect structures we find lack 
of agreement between models and experiments. Experi- 
mentally the study of stripe formation has taken a sub- 
stantial step forward with the work of Harrison et al [10] 
on carefully prepared two dimensional diblock copoly- 
mer systems which order into striped systems appearing 
to be two dimensional smectics. For quenches into the 
appropriate temperature range they find ordering which 
proceeds, in the scaling regime, via the process of the 
annihilation of a set of disclination quadrapoles. They 
also find that characteristic growth laws grow in time 
Lit) w t x with exponent x = 1/4. 

The Swift-Hohenburg (SH) model (see Eq. (11) be- 
low) [11] represents a simple model for producing stripe 
pattern growth. In a previous paper [8] we discussed 
the distribution of defects (grain boundaries and point 
defects) which govern the ordering in the SH model in 
two dimensions. We found [8] the kinetics dominated by 
grain boundaries with a small number of free dislocations 
and a smaller number of free disclinations. Thus the sim- 
plest model, the SH model, does not produce the defect 
structure seen in the experiments. 

In this paper we present an alternative model de- 



scription, the nonlinear phase (NLP) model, for grow- 
ing stripes. It is motivated as an approximation to SH 
model. However the model grows stripes via a quite dif- 
ferent defect structure compared to the SH model, but 
shares with it the characteristic feature of grain bound- 
aries. This model could also, as with the SH model, be 
constructed by appeal to symmetry, simplicity and ana- 
lyticity. 

There appears to be a variety of pathways to equlib- 
rium in these systems. One may need to appeal to more 
than just symmetry and analyticity to obtain a quanti- 
tative description of the ordering process in striped sys- 
tems. 



II. THE NONLINEAR PHASE MODEL 

The nonlinear phase model can be obtained as an ap- 
proximate phase-field model for the SH model defined 
by 



(1) 



where ip is a real scalar field. We only consider the zero 
temperature case in this paper, so there is no noise term 
in the above equation. Assume that the solution of the 
SH model can be written in the single mode form [12]: 

V>(x, t) = ip cos(k(x, t)) (2) 

where the amplitude ipo saturates quickly at the ground 
state amplitude \/4e/3. Substituting this ansatz into the 
SH equation, ignoring amplitude fluctuations and higher 
harmonics, we find that the coefficient of ij)Q sin(«;(x, t)) 
satisfies 



where 



-k = Vj [V 2 Qj + 2(1 - Q 2 )Q ] 



Q(x,<) = V/e(x,i) 



(3) 



(4) 



is the local wavenumber of the stripes. Eq. (3) can be 
written as 
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If we take the gradient with respect to Vj, Eq. (3) can 
be written in the form 



dtQi = V«V 



SH E 

8Q, 



(7) 



where the driving free energy is of the standard Ginzburg- 
Landau form 



n E [Q] = Jd 2 x 

= J d 2 x e(x) . 



^(VQ) 2 + i(Q 2 -l) 2 



(8) 



Except for the fact the Q is longitudinal and one has an 
anisotropic diffusion tensor (VjVj rather than V <J»j), 
this would just be the TDGL model for a conserved vec- 
tor order parameter. The model (nonlinear phase model 
or NLP model) defined by Eqs. (2), (3) and (4) corre- 
spond, as we now show, to a well defined stand alone 
model for producing striped patterns. While we have 
roughly derived this model from the SH model, its or- 
dering defect structures are quite different from the SH 
model. 



III. NUMERICAL RESULTS 
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FIG. 1. cosk > at t = 1000 in a 512 x 512 system. 

This is a fairly early time for this system and we see that 
we have generated a set of local donuts. We have many 
concentric circles with few direct paths through the sys- 
tem. We clearly have layers but they are strongly bent. 



A. Qualitative Results 

Let us begin with a qualitative description of the pat- 
terns generated by this model. We put Eq. (3) on a 
lattice with spacing Ar = 0.7854 and time step 0.01. We 
use the isotropic form for the Laplacian given by 



M) 



(Ar) 2 



2^1^ 10 
3^ + 6 ^ 3~ 

NN NNN 



(9) 



where NN and NNN mean the nearest neighbors and 
next-nearest neighbors respectively. <p can be the order 
parameter n or the component of Q. (i, j) is the lattice- 
site position. We run these equations on lattice of vari- 
ous sizes using initial conditions where the phase variable 
is chosen to have a value randomly distributed be- 
tween -k and k where k = 0.01/V2 = 0.0070711. 

In Fig. 1 we plot those points where cos k is positive 
for a system, which is on a 512 x 512 grid, after an evolu- 
tion to time t = 1000. The stripe wavenumber Q 2 orders 
very rapidly so that by the time t = 1000 we have an in- 
terconnected pattern of basically compact objects. These 
appear to be nucleated objects with radial symmetry. 




FIG. 2. cosk > at t = 5000 in a 512 x 512 system. 
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When we reach a time of 5000 we see that most of the 
donuts have opened and stripes are forming and winding 
through the system. Clearly we see that the remaining 
centers of the donuts are serving as cores for large discli- 
nations. 
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FIG. 3. cos k > at t = 15000 in a 512 x 512 system. 

This pattern coarsens as one moves to t = 15000, as 
shown in Fig. 3, 
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FIG. 4. cos k > at t = 25000 in a 512 x 512 system. 



and Fig. 4 where t — 25000. In Fig. 4 one sees a target 
pattern in the upper left. One expects this to eventually 
break open. 



These patterns can, at later times, be rather com- 
pletely characterized by their defect structure. The de- 
fects can be found by looking for positions where the am- 
plitude Q 2 is significantly smaller than its ordered value 

Ql = i- 
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FIG. 5. cosk > at t = 15000 in a 512 x 512 system. 
Those points with a small Q 2 are also shown. 



In Fig. 5 we plot those points where Q 2 < 0.5 for the 
pattern in Fig. 3. We see that we have a rather complete 
map of the defect structure seen in the layer pattern. The 
small amplitude lines (SAL) network defines the pattern. 
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FIG. 7. The Q field for the pattern in Fig. 6. The vector 
on each site is shown. 



The +1 vortices are compact but the —1 vortices are 
not. This can be seen in Fig. 7, where some of the — 1 
vortices occupy entire edges. So the defects in the system 
are more like a domain-wall network rather than a set of 
vortices. The domain walls correspond to the edges in 
Fig. 6. 



FIG. 6. The pattern cos k and k(x) for a 64 x 64 system at 
t = 800. 



B. Quantitative Results 



In Fig. 6 we plot the pattern cosk > and k(x) for 
a 64 x 64 system at t = 800. Clearly the peaks and 
anti-peaks correspond to the target centers of the pattern 
cos k 7 and the edges correspond to the SAL network. We 
also observe that the peaks and anti-peaks correspond to 
the +1 vortices in the Q field, and on some of the edges 
there is a —1 vortex of the Q field, as is shown in Fig. 
7. The numbers of +1 and —1 vortices are equal, while 
the numbers of the peaks (anti-peaks) and the edges are 
usually not equal, so not every edge has a —1 vortex on 
it. 



At long times after a quench, phase ordering systems 
typically enter a scaling regime with a single character- 
istic length L. Here we check scaling for the correlation 
function C K (r, t) = (k(x + r, t) k(x)}. In Fig. 8, we 
show this correlation function at 8 different times. The 
data is scaled so that we can see whether it obeys a scal- 
ing law. The correlation length can be extracted from 
C K (r ,t)/C K (0,t) = 1/2 where r cx L. The result is 
shown in Fig. 9. Since |Vk| ^constant, we expect that 
k <~ L and C K (0, t) — (k 2 ) ~ L 2 . Thus we expect the 
scaling form C K (r,t) = L 2 F{r / L). The direct measure- 
ment of (k 2 ) is shown in Fig. 10. From Fig. 9, we get 
L ~ t 60 . From Fig. 10, we also obtain L - i 60 . 
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FIG. 10. (k 2 ) oc t 1 - 20 . So L 2 ~ t 1 - 20 
aged over 41 runs. 



and L ~ t ' 60 . Aver 



FIG. 8. The correlation function C K {r,t) for the 256 x 256 
system is plotted, for 8 different times, in scaling form. The 
distance r is scaled by L Ha lf = fo given in Fig. 9. At least at 
small r, the function has a scaling form. The data is averaged 
over 41 runs. 
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We measure next the number of defects in the system as 
a function of time after the quench. In Fig. 11 we plot 
the density of sites where Q 2 < 0.4 in the scaling regime. 
This curve is well fit by 



n v = at 



(10) 



with a = 3.89, n = 0.506. Since the total number of the 
defects is proportional to the area of domain walls, the 
number density is proportional to w ■ L/L 2 = w ■ L~ x , 
where w is the average width of the domain walls. Thus 
we estimate L/w <~ 0.5. 




FIG. 11. Density of the defect points where Q 2 < 0.4. Av- 
eraged over 38 runs. 



FIG. 9. We extract the correlation length from 
C K (ro,t)/C K (0,t) = 1/2. ro = Lnaif is proportional to L. 
We find L ~ t 0,60 . Averaged over 41 runs. 



It will turn out to be necessary to allow the width w to 
be a function of time. 

We need an independent method for determining L or 
w. We do this by employing another method to measure 
the length of the defect network directly. Basically this 
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is a coarse-graining method. We put the defect network 
on a lattice with lattice spacing 8Ar and count the num- 
ber of sites it occupies. This number reflects the length 
of the network much better than the previous measure- 
ment. Other sizes of the lattice spacing can also be used 
as long as the grid size is larger than the width of the 
domain walls and much smaller than their lengths. The 
length density of the network is shown below and can be 
fit to <~ i~ - 60 . Since the length density is proportional 
to L/L 2 = L^ 1 , we get L <~ t 0,60 , which is consistent 
with the measurement of the correlation function. 
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FIG. 12. Length density of the defect points where the 
length of the network is in fact the number of sites that the 
network occupies in the lattice with larger grid. The dashed 
line is proportional to t~ 60 . Averaged over 38 runs. 



We can then conclude that the average width of the 
domain walls w ~ t 010 . In the time regime we study, the 
width w is about 2^4. We don't think w will increase 
for ever, it may stop increasing at some later stage. But 
before that can happen in our system, finite size effect 
enters. 

The growth laws for L and w can explain the growth 
exponents of other quantities. We give two examples be- 
low: {Q 2 } and the energy density e in Eq. (8). 

The ordering of the field Q is characterized by the av- 
erage over all sites of Q 2 (z,j). We obtain the results 
shown in Fig. 13. It can be seen in the figure that there 
are two regimes where the data can be fit to a form 



(Q 2 ) = 



bt- 



(ii) 
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FIG. 13. (Q 2 ) v.s. t. Averaged over 20 runs. 

In the shorter time regime 1000 < t < 7000 the data 
is fit with with a = 0.930, b = 3.15, and n = 0.45. In 
the longer time regime 7500 < t < 10000 we have the 
fit a = 0.950, b = 2.40, and n = 0.511. In both fits the 
exponent n is near 1/2. We also see that in both fits 
the final value of Q 2 is larger than 1. Clearly more work 
needs to be done to establish the nature of this cross over. 
It seems likely that the final asymptotic value of Q 2 is 
greater than one due to finite size effects and the freezing 
of grain boundaries for long times. 

In this system, Q 2 is approximately 1 away from the 
domain-wall network and small on the domain walls with 
an average c < 1. Because the average area s of the do- 
main walls in one domain is proportional to wL, and the 
domain's area is A ~ L 2 , we have 



(Q 2 ) 



> 



s) ■ 1 + s ■ c] = 1 - (1 



"•3 



(12) 



where s/A cx wL 1 . So the average of Q 2 has the form 
of (Q 2 ) -Ql~ w/L. Thus we can identify L/w - i 5 . 
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FIG. 14. The energy density of the system. Averaged over 
20 runs. 

In Fig. 14 we show the path toward equilibration of 
the effective energy density e defined by Eq.(8). This re- 
sult seems to be in agreement with that for Q 2 . We have 
a good fit to 

e = a + bt' n (13) 

with a = 0.0056, b = 2.526 and n = 0.5288. Since the 
energy above the ground state should be proportional to 
the area of the domain walls, the energy density is pro- 
portional to wL/L 2 — wL^ 1 . Again we get L/w ~ t - 5 . 



ponents [8] . Our guess is that in SH model there are also 
more than one characteristic length. 

In conclusion we see that models with the same symme- 
tries can vary significantly in the defect structures pro- 
duced during ordering. The search continues for mod- 
els where one can dial the relative abundance of grain 
boundaries and free dislocations and disclinations. The 
ultimate goal is to match the models with a given exper- 
imental system. 
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IV. DISCUSSION 



In the NLP model we know that there are analytic 
vortex solutions related to those for the XY model. This 
would seem to favor coarsening via a set of isolated vor- 
tices which pair up and annihilate. Instead we find large 
— 1 vortices forming a domain wall network. 

This system does not appear to generate dislocations 
and, as such, is quite different from the SH model which 
has a significant density of dislocations. The nonlinear 
phase model, on a larger scale, is growing targets and 
disclinations. 

The NLP model introduced here helps to deepen our 
feeling that we do not have a good understanding of 
the general mechanism of stripe formation. In the SH 
model and the appropriate experiments the ordering is 
slow compared with the simple model of point defect an- 
nihilation. In this NLP model the ordering is "faster" 
than the simple model. 

In this model we find two characteristic lengths, L and 
w, which together explain the different exponents we ob- 
served. In the SH model, we also observed different ex- 
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